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Abstract 

We reexamine the solution of the minimal Zee model by comparing with 
the data of the SNO experiment, and conclude that the model is strongly 
disfavored but not yet excluded by the observations. Two extensions of the 
Zee model are briefly discussed both of which introduce additional freedom 
and can accommodate the data. 
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I. INTRODUCTION 



To understand why neutrino masses are non-zero is one of the most important subjects 
in particle physics. Assuming only left-handed neutrinos UaL (a = e, /i, r) the minimum 
standard model predicts m{va) = since adding a bare Majorana mass term Mab^aL^bL 
with Mab = Mba for these neutrinos violates gauge invariance (because of the concomitant 
^aL^bL - type terms) and hence renormalizability. Thus, further states must be added to 
accommodate the neutrino masses, and we look for the greatest economy and simplicity in 
doing this. For a review of the theory, see 

II. MINIMAL ZEE MODEL 

The Zee model is one of the most economical possible scenarios. In the Zee model, 
the Majorana neutrino masses are generated by a one loop diagram. The origin of the 
smallness of the masses come from this feature. Hence one of the present authors has 
discussed the comparison with the neutrino experimental data in Ref. [^. And several 
other authors have also discussed it and shown similar results [^Q. The solution of the 
neutrino mixing is bimaximal { 6i = 63 = 11/4,62 = with 6i as defined below). This 
agrees with the atmospheric neutrino data. For the global analysis of solar neutrino data 
0, this solution corresponds to the large-angle MSW (LMA) solution or the just-so vacuum 
oscillations (VAC). However the result from the Zee model does not agree so well with the 
LMA solution from the recent analysis |^ and the recent data of the super-Kamiokande 
disfavors the VAC solution. So we need to reconsider the possibility of the LMA solution in 
the Zee model. 

Recently the SNO group announced their experimental data of solar neutrino fluxes 
from decay measured by the charged current reaction rate, which is 0'"'"(z/e) = 1-75 ± 
0.071 o;ii ±0.05 X lO^cm ^. By combining this with the data from the Super-Kamiokande 
and using the values of the total flux expected in the Standard Solar Model (SSM) the 
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survival probability of z/g was reported in ref. 0. It is 

A.CC 



JSSM 



0.347 ± 0.029+°:°^^ 



(1) 



where the first error is from SNO and the second from the SSM(Solar Standard model) theo- 
retical error. Using a combination of (tf^{i'e) and (j)^^{vx) from SNO and Super-Kamiokande 
leads to the estimate IPI 



0.322 ±0.076. 



(2) 



This suggests the survival probability P{ve — >■ Ve) is nearer to 1/3 than 1/2 and this is new 
information for us to analyze. 

In this paper, we rediscuss the compatibility of the answer from Zee model with the LMA 
solution [0 by comparing with the recent SNO data. 

Under the Zee anzatz, the neutrino mass matrix in the flavor basis (e, /i, r) is 



M 



^ M,^ Mer ^ 



M,^ M^, 



Me, M.r 



^ mi ^ 



U 



J 



V 



ma 
ma 



(3) 



where mi,m2,m3 are the eigenvalues of Ai and U is the unitary matrix to diagonalize it. 
Ai is real, traceless and symmetric. From the tracelessness condition, 



mi + m2 + ms = 0. 



(4) 



This condition is a strong constraint. The mass pattern of exact solutions which satisfy the 
atmospheric neutrino data is 

mi = -m2, ms = 0, (5) 
and the ratio between the two neutrino squared-mass differences is r = A^/A^ = \mf — 



m2|/ |m| — m; 



0, where subscripts s, a refer to solar, atmospheric respectively. We will 



examine r > later. 
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With this situation, the allowed mixing matrix is the bimaximal one with 6i = 11/4,62 = 
and ^3 = 7r/4, where the definition of the mixing angle i^. 



U 



\ 



C2C3 C2S3 S2 

-C1S3 - S1S2C3 C1C3 - S1S2S3 S1C2 
S1S3 - C1S2C3 -S1C3 - C1S2S3 C1C2 

with Si and q standing for sines and cosines of 6i and the bimaximal mixing matrix is 



V 



(6) 



/ 1 
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1 _ 1 J_ 

V 2 2 / 



(7) 



To discuss the neutrino fiux from the sun, we have to solve the neutrino propagation equation 
in the matter as follows: 



dt 



( \ 
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( \ 
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where A = 2\/2GpN(,E, N^, is the density of electron neutrino in the sun, E is the energy 
of the neutrino. On the condition of eq.(5), |mi| = \m2\ and = 0, 



ml + A 







V 








12 12 

1212 



(9) 



Then, the rotation among the weak eigenstates and the mass eigenstates in the center of the 
Sun z/^, u^p) can be expressed as follows: 



V^3 / 



1 
^ ^ 

V~72 ^ 73/ 



^ . .rn ^ 



V^3 / 



(10) 



^There is no CP violation in Zee model. 
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where we have taken the hmit of large electron neutrino density, {A oo). At t = an 
electron neutrino is produced in the sun and it is composed mainly of the state z/^ in the 
hierarchy ^ we are considering in this work. 

ke(0) >= > . (11) 

The time evolution of this state to time t is 

\iy,{t) >= e'^o^'^'luip > (12) 

where A2 is the eigenvalue of for i/™ state. Since the neutrino is measured on the earth, 
the state of z/g is expressed by the mixing Eq.(|^). The amphtude is 

<Z/e|z/e(t)>=^e*^-'*. (13) 

Hence the survival probability is 

P(z/e ^ Z/e) = I < Ue\Ue{t) > T = ^- (14) 

This is the result from the exact Zee anzatz with r = and is significantly disfavored by 
SNO data. We are led to consider a more realistic case, r 7^ and \mi\ 7^ |m2|,m3 7^ 0. 
For the more realistic case, we rewrite the mixing matrix of Eq.(|^) by using the following 
parameters which show the discrepancy from exact bimaximal mixing. 



Cl = COs(--^l)~-^(l +£!--£?), (15) 

vr 1 1 

si = sin(- -^1) ~ -^(1 -£i - -£?), (16) 

C2 = cos(6) ~ (1 - lei), (17) 

S2 = sin(^2) ~ £2, (18) 



^As we show later, in the case ms 7^ 0, the mass hierarchy we need is m2 > mi ^ 7713. By the 
matter effect, the hierarchy in dense matter is m2 ^ nii ^ 771,3. So we can define that the electron 
neutrino at the production point is composed mainly of the state 1^2- 
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,7r 



C3 



S3 



cos(- -^3) 
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sm 
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1 . 1 2x 

^(1+^3-2-3), 

/3J, 



£^3 



(19) 
(20) 
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where we neglected e"^ and is the difference of the angle from the bimaximal case and 
Ex = sin^j,'. By using this expansion, we can find the following relations up to 0(£:^) among 
the parameters from the conditions that the diagonal element of Ai are zero, 



7711 + 1712 = -263(7171 - 7712) 



62 — oei^s. 



(21) 
(22) 



By these relation and the tracelessness condition mi + m2 + ^3 = 0, we find the relation 
between r = A^/Aq and £3 to be 

863 



12 21 
\7n\ — 

1 2 2T 
\m\ — 7n^\ 



l-4£3-12£i" 



(23) 



The behavior of r versus £3 is shown in Fig.l. The experimental data suggest that the ratio 
r satisfies r < 0.1 at 90% confidence level. So to satisfy this upper bound on r, £3 cannot 
take too large values. From Fig.l, we find the magnitude of £3 is smaller than about 0.015. 

From the neutrino propagation equation, Eq.(^), after replacing by mg = — 2e3A, 
which comes from Eq.(21) and the traceless condition, we find 



where 



(24) 



Ml 



\ 








4^ 



4^ 



4^ 4^ 



(25) 



and 



Ml 



-ir2^A + e,A + e2T.) 
\ ^(A-^iA-.^S) 



-i72(^ + ^i^ + ^2S) ^(A-.iA-.^S) 

i(£lS-3£3A) -i£3A 



i(£iS-3£3A) 



(26) 



in which S = mf + m"^ and A 



m|. In the case of £3 7^ and £1 = £2 = 0, the 



mixing matrix which diagonahzes in the matter is 



where 



/ 



V2 + K^ 



V2 K 



— ^ 1 

-1 



V2 



V2 ) 



-A - 2£3A - MA + 2£3A)2 + A^ 



(27) 



(28) 



This means that the mixing angle ^™ satisfies the following relation, 



cos 9\ 



sin 



(29) 



^The atmospheric neutrino data favored ei = and from Eq.(22) £2 = is also favored. 
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FIGURES 
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FIG. 1. r as a function of £3. 




FIG. 2. The survival probability of the electron neutrino. The horizontal lines show the bounds 
from the SNO experiment: the solid lines exclude the SSM theoretical error while the dotted lines 
include it. The diagonal line is from Eq. pO|) with Pc set to zero. 
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The production point of the neutrino from decay is in the central core of the Sun so that 
the electron density A has a very large value. At the production point, sin0™ ~ 1. This 
means the electron neutrino at the production point is composed mainly of 1^2 as for the 
earlier discussion on the case ms = 0. Then the survival probability is 

< P(z/e ^ Ue) >= ^ - ^3(1 - Pc) (30) 

where Pc is the hopping probability at the resonance point and almost zero in the LMA 
region. The dependence of £3 on the survival probability is shown in Fig.2. The constraint 
from SNO experimental data is given in Eqs.(|I]) and (0). 

To satisfy this condition, the parameter, £3, would need to be larger than 0.1. But then 
the value of r becomes significantly greater than 1 and this is very hard to accommodate 
within the data on and A^. The data suggest A^ > 1.5 x lO^^e^^ |[Tl[] while A5 < 
2 X IQ-^eV^ ^MM at 90% C.L. 

We conclude that the minimal Zee model is strongly disfavored but not yet fully excluded 
by the SNO/SuperKamiokande results. For the minimal Zee model to survive, one would 
need A 5 to get larger than A^. 



III. EXTENSIONS OF THE MINIMAL ZEE MODEL 

Before looking at specific extended Zee models we consider adding diagonal elements in 
the neutrino mass matrix and relaxing the tracelessness condition to the more general 

mi + 7712 + rris = 5 x mi (31) 

where S is the shift from and normalized by mi. Then, we request the Mee = in order 
to avoid the difficulty from double (3 decay. We then find the constraints as follows: 

mi + m2 = -263 (mi - m2), (32) 
4^3 + 5 - 2^35 

nis = mi. (33) 

i — ZEs 
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From these conditions, 



rrii — 1712 



nil ~ ""^3 



863 



{1 - 6^ - 4e,{l - 6y - 4el{3 - 6){1 ~ 6)} 



1 - 2e;V 



(34) 
(35) 



To satisfy the survival probabihty result requires that £3 is greater than about 0.1, while r 
must be smaller than 0.1. However this condition is not satisfied unless we keep the mass 
relation mi > m^. The correlation between the ratio r and m^/mi when £3 = 0.1, 0.15 
and 0.2 is shown in Fig. 3. In this case, to satisfy the conditions from experiments, the mass 
hierarchy is |m3| ~ \Smi\ = 2|M^^| > |4mi| > \m2\ > and the parameter 6 has to be 
larger than about 4 from Fig. 3. So a very large deviation from the tracelessness condition 
is needed. 




FIG. 3. The correlation between r and ms 

We consider the case that the all diagonal elements of M are not zero. Here we param- 
eterize the diagonal elements as follows. 



Mee = See X TT^i, 



(36) 
(37) 



Then, by the angles we defined in Eqs.(|T5D - (^OD, the relations Eqs.(0) and (|3^ are changed 
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to 



From these relations, 



And then, 



1711+7712+7713 = {6ee + 25^^)mi, 

mi + 7712 + 2£:3(mi — 7712) = 26ee X mi. 
ms = 2e3(mi - 7712) - {See - 25^^)mi. 

4(l-4e)(2£3-4e) 



(38) 
(39) 

(40) 

■ (1 - 2^3)2 - {(1 - 2£3)(5ee + 25^^) - 2(2£3 " 6ee)V' ^^^^ 

In Fig. 4, A5 = rAa, with = 3 x lO^^eV^ is plotted versus 6ee using Eq.(|?T|). To reduce 
r, 5ee should be near to 2£:3( or to 1), because the numerator of eg. (PI) has the factor of 
(2£:3 — 6ee)- Heuce, this condition will admit consistency with the LMA solution. However 
we have to note that 6ee is constrained by neutrinoless double P decay experiment. In this 
case (5^^ = 0), the mass pattern is 

mi ~ —7712, nT'3 ~ Mee ~ 2e3mi, (42) 

and mi ~ O{10~^)eV. This is not yet excluded by neutrinoless double f3 decay experiment. 

In this discussion, we found the element Mee is important to realize the LMA solution 
from the model based on the Zee model. 

Ag = 3 X IQ-'^eV^ 
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FIG. 4. As plotted versus 6ee using Eq.^ for a fixed = 3 x 10'^ eV'^. 

In the light of these remarks we now consider briefly two extensions which introduce more 
freedom and can accommodate all the present data: (1) More than one independent Higgs 
doublet coupling to charged leptons, enabling off-diagonal flavor vertices for those doublets 
without a vacuum value (2) Addition of a singlet doubly-charged scalar ||T^|-0. 

The neutrino Majorana mass in extended model (1) arises from the graph in Fig. 5. 



/i < V, > 




FIG. 5. Majorana neutrino mass for the minimal Zee model or for the extended model (1) 



7 



If 
I 

I 

I 
I 

I 



V 



FIG. 6. /i ^ 67 in extended model (1) 
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The Yukawa couplings in this model are 



(43) 



where is left-handed lepton doublet and Ir is right-handed charged lepton. The 0j are 
Higgs doublets. After SU(2) symmetry breaking the remaining charged Higgs coupling is 



{sindHt - cos9H+ 



(44) 



where 6 is mixing angle between the charged Higgs and the singlet Zee scalar, m* is the 
mass of charged lepton, Vi are the vacuum expectation values for the Higgs bosons and 
are mass eigenstates. The g2 has flavor off-diagonal elements and by this feature diagonal 
neutrino masses arise from Fig. 5. 

The off-diagonal couplings also contain a /i — > 67 flavor violating decay [1^ arising from 
the diagram of Fig. 6. The branching ratio is 

a 



BR{n 67) oc 



vi 



+ 



-\ 2 



Ml 



Ml 



(45) 



where g2 is the Yukawa coupling with flavor off-diagonal elements. The experimental bound 
is BR{j2 67) < 1.2 X 10"^^. From this bound. 



M4 



< 2 X 10"''(G'e1^) 



(46) 



where ^ 



'""t^ + If we take M ~ lOOGeV, this imphes 



c/2 < (4 ~ 7) X 10 



-3 



(47) 



While, the mass elements are 



j2 _j_ y2 

Seetni = OC — ^/ifi'2/eM"^/., 



nil 



_ r 2 



Vl 



r- - OC llfeuVn,,- 

^/2 ^ Vl ^ 



(48) 
(49) 



^We neglected the contribution from the first term in Eq.(44) in this discussion because it is 
already smaller than the §2 term if g2 ~ 0(10^'^). 
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where /e^ is the Zee scalar couphng and we neglected the term of fert^r because /er/ fev oc 
m^^/rn^ from |Me^| ~ |Me^|. The ratio is 

v^5ee ~ ^ OC ~ "^g,. (50) 

We need this ratio to be near to 2v/2e3. So we find: 

~ 2V2em,-^^ ~ 0(10-^). (51) 

Ui -t- 

estimated for the case Vi = v^- This condition will then comply with the bound from /i — > 67. 
Then extended model (1) is consistent with all the data. 

Finally we consider extended model (2) which contains new couplings with a doubly- 
charged scalar singlet k^^: 

hxxllCLlk-^+ + h.c. + K{h+h+k— + h-h-k++), (52) 

where is a doubly charged SU{2) singlet scalar which couples to only right handed 
charged leptons. The contribution to the Majorana neutrino mass now comes from the two 
loop diagram of Fig. 7. 



/ Ay , \ 

' mi I mi 
> 1 — < >< * — \ — » — * — I 4 

h Ir Ir II 

FIG. 7. Majorana neutrino mass from two loop diagram in the extended Model (2). 



The diagonal mass element is 



^^i/i— f/eu5 (53) 



14 



where mx {X — /i, r) is the charged lepton mass, is the mass of doubly charged scalar 
and F show some log function. Here we neglected the term /^^ for the same reason we 
neglected it in Eq.(48). To get the required mass element, Mee has to be around 2e3mi. 
Namely, the ratio between the diagonal and off-diagonal elements of the neutrino mass term 
should be near to 263. 

We may always realize this condition which depends on unknown parameters, and hence the 
extended model (2) is consistent with all the neutrino data. 



IV. DISCUSSION 

The minimal Zee model is very economical as a simple way to introduce neutrino mass 
into the standard model. However, it is seen to be barely consistent with the combination of 
SuperKamiokande and SNO data. It would need A5 > which looks like a considerable 
stretch from the observations at hand. 

On the other hand, if we enrich the model by either (1) adding further Higgs doublets 
couphng to the charged leptons and to the singlet charged scalar or (2) by adding a doubly- 
charged singlet scalar, there is enough freedom to accommodate the SNO data. 
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